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To  be  realistic,  any  mode!  of  structural  fatigue  should  allow  for  atlrlilon  due  to 
collapse  andjor  war  damage,  hijacking,  etc.  In  OHe~crack  structures,  attrlikm  adds  a slugle 
term  to  the  Fokker-Planck  equation  which  statistically  describes  cruek  growth. 
This  remarkably  simplifies  the  reliability  computation,  and  elucidates  the  interoetkm  of 
fatigue,  inspection,  ultimate  loads  and  other  risks. 
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Cnck  size 

Crack  ii»  at  initiation  time 
Critical  crack  size,  It(e)  infinite 
Expectation 

Conditional  probability  density  of  crack  length  (see  ^.)> 

Probability  distribution  of  crack  length 

Density  of  initial  life  t or  t. 

Distribution  of  initial  life  t 

Attrition  or  hijack  factor  of  initial  life  survivorslup  fiinction 
Generic  form  of  moment  generating  function,  £.exp(u.) 
Crack  growth  rate  (local  average) 

Risk  fiinction 

Risk  in  the  absence  of  a crack  ‘‘hijack  risk” 

Equivalent  time  or  cycles 

Initiation  time 

Elapsed  time  since  initiation 

Dominating  iime  used  in  analysis  of  runaway  cracks 

Marginal  density  and  distribution  of  crack  length  at  failure 

Conditional  density  of  crack  length  at  a given  life 

Arithmetic  mean  life  to  initial  failure 

Symbol  for  set  intersection,  used  here  as  “and” 


1.  INTRODUCTION 

Since  the  early  1960s  research  in  fatigue  prediction  has  divided  into  the  distinct  areas  of 
crack  growth,  reliability  or  inspection  and  finally  structural  fatigue.*'*'*  Until  now,  most  of  the 
models  assumed  for  reliability  have  either  been  forced  to  make  many  simplifications*  which 
have  been  criticised*'*  or  else  the  calculations  have  been  inordinately  lengthy  and  possibly  ill- 
conditioned.  These  remarks  apply  more  particularly  to  the  assumptions  used  when  inspections 
are  included.  Apart  from  this,  many  models  involve  structures  or  data  that  are  highly 
idealised.*'*'* 

Structural  fatigue  is  now  seen  to  involve  the  statistical  interaction  between  initiation  (damage) 
and  crack  growth  at  several  fatigue  critical  regions  in  a structure.  This  concept  developed 
naturally  from  the  separation  of  fatigue  into  distinct  phases  of  damage,  cracking  and  final  failure 
which  is  now  the  subject  of  reliability  theory.  The  buic  problem  is  now  to  find  the  multivariate 
distribution  of  crack  lengths  at  any  time.  Although  the  general  theory  is  multivariate,  the  single 
crack  case  is  not  trivial  and  forms  usefUl  illustrations. 

In  the  earlier  thinking  of  Ford,*  the  crack  length  distribution  was  used  to  find  risk  rates  for 
reliability  estimates  in  a separate  but  simultaneous  computation.  Thus  until  now  general  theory 
of  structural  fatigue  has  made  no  allowance  for  failures  or  retirements  after  inspection.  In  this 
paper,  such  an  allowance  is  made  for  sin^e  cracks.  In  terms  of  structural  fatigue  theory  the 
extension  is  almost  trivial  but  the  reliability  theory  is  profoundly  simplified  and  generalised.  The 
success  of  this  one-crack  model  means  that  prediction  of  structural  fatigue  must  now  intinuitely 
involve  reliability  or  attrition,  so  that  a multi-crack  reliability  theory  is  now  needed.* 


2.  GENESIS  OF  NEW  ONE-CRACX  MODEL 

The  long  fatigue  life  normally  expected  lends  itself  to  measurement  in  terms  of  pseudo- 
times such  as  flights,  service  hours,  years,  etc.  These  are  pseudo-times  because  they  indicate 
length  of  service  rather  than  physical  time.  Since  practical  time  units  contain  many  fatigue 
cycles  it  is  mathematically  convenient  to  let  time  be  continuous,  the  common*'*'*  but  implied 
practice  which  we  also  follow.  With  such  an  idealisation,  some  events  may  be  slightly  anomolous 
in  terms  of  discrete  cycles. 

We  make  a distinction  between  retirements  after  inspection  at  fixed  times  (termed 
“inspection”  for  brevity)  and  “attrition”  at  arbitrary  times  from  a continuous  combined  risk  of 
overload  failure,  war  damage,  hijacking,  etc. 

Following  previous  definitions,*'*'*  fatigue  is  now  assumed  to  proceed  in  three  distinct 
stages.  The  first  of  these  ends  with  the  initiation  (at  length  a,)  of  a crack  whose  subsequent  length 
a is  aMumed  herein  to  be  a deterministic  function  of  elapsed  time  in  the  second  stage.  The  third 
stage  is  attrition  (including  overload  failure),  and  the  structural  fatigue  problem  is  to  find  the 
corresponding  probability  density  function  ^0  of  life,  which  may  allow  for  inspection.  Because 
initiation  is  determined  by  cumulative  damage,  the  pre^rack  stage  may  also  be  termed  the 
damage  stage;  “damage”  does  not  refer  to  cracks,  except  that  the  next  stage  begins  with  a crack 
(of  length  a»). 

These  stages  are  described  by  the  three  given  functions /(0<  R(n)  and  rfa)  relating  to  initiation, 
crack  rate  and  risk.  In  practice  of  course,  they  are  the  output  of  preliminary  estimations  using 
relevant  data  and  any  fatigue  method  whidi  averages  the  effects  of  load  sequence.  If  inspections 
are  included,  an  operating  characteristic  F(a)  is  also  required  relating  success  of  inspection  to 
defect  size  (see  page  9). 

Thus,  the  model  developed  below  produces  a life  distribution  4(0  from  the  data  fiinctions 
/(O,  R(n)>  >(a)  and  possibly  ^a).  It  should  also  be  noted  that,  like  /,  the  crack  length  a may  be 
arbitrarily  defined  as  long  as  it  adequately  describes  the  cracked  structure.  For  example,  crack 
length  or  area  may  be  used,  but  the  equt%  valid  log  (area)  may  be  more  convenient. 


Crick  Imgth  • 
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If  there  U only  one  crack  the  damage  is  unaffected  by  crack  growth  for  a given  structure  so 
that  initial  failure  can  be  immediately  predicted  from  an  appropriate  damage  rule  or,  if  inspec- 
tions are  of  main  interest,  experimentally  observed  intiations  may  be  used.  If  there  is  no  inspec- 
tion and  no  final  failure  then,  with  deterministic  cracking,  any  crack  that  begins  between  any 
two  initiation  times  retains  the  same  probability  measure  from  those  times  onwards  (a  runaway 
crack  is  a special  case  of  this). 


Inspection 
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Density  at  t cycles 
/(ajt) 
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FIG.  1 EFFECT  OF  INSPECTION  ON  CRACK  LENGTH  DENSITY 


In  Figure  1 however,  if  there  is  an  inspection  at  t the  probability  of  having  the  crack  between 
A and  B will  be  discontinuously  reduced  at  that  time  by  an  amount  that  depends  on  oa,  an  the 
poMible  lengths  at  that  time.  The  overall  rejection  rate  is  an  average  of  such  quantities. 

Although  deterministic  cracking  is  assumed,  local  randomnem  can  be  included  in  a quasi- 
deterministic  way  and  this  has  been  done  in  the  extensions.  In  the  next  section  the  possible 
cracks  are  also  autonomous,  i.e. 


db/dk  - Kfl) 


regardless  of  the  time.  If  this  is  integrated  then  the  inversion  of 


f'A.,. 


»e-  U 


indicates  that  any  crack  length  in  this  one  dimensional  case  depends  only  on  the  growth  time  u. 
Fignn  2 shows  some  typical  ilhistrationt  of  this.  However,  alt^^  this  is  important  in  general 
structural  fatigue  it  is  not  necessary  here.  The  more  gens^  case  is  treated  as  an  extension  so 
that  the  first  presentation  is  not  too  chitteied. 


L2  AttiMw 


If  iiiq)ectioiu  are  included  in  this  manner,  it  is  only  fair  that  the  continuous  unscheduled 
“inspection”  by  the  environment  is  also  included.  That  is,  the  probability  of  overloads  leading 
to  failure  should  be  covered;  this  is  the  subject  of  conventional  reliability.  As  used  below  it  will 
include  the  continuous  risk  of  removal  from  any  cause  of  interest,  and  the  risk  on  an  uncracked 
structure  r«  will  assume  a particular  importance.  The  latter  will  affect  the  probability  density  of 
initiation  so  that  thus  must  also  be  included  in  the  model. 

In  order  to  simplify  the  language  it  is  convenient  to  introduce  the  term  attrition,  or 
(generalised)  failure,  for  the  process  described  by  reliability  theory.  Thus 

Risk  function  - Conditional  probability  of  attrition. 


2,3  Notational  Caavcatioas 


In  the  following  a generic  notation  will  be  used  so  that /(.)  and  F(.)  will  refer  to  the  (wobability 
density  and  distribution  functions  of  the  arguments  and  the  form  ot  the  function  will  differ 
with  the  argument.  As  a special  case  ^.)  refers  to  any  density  affected  by  attrition,  including 
life  distribution.  Af.(u)  will  refer  to  the  moment  generating  function  of  the  random  variable 
indicated  by  this  subscript. 

For  example. 


M/u)  = 


e“*  dFW. 


(0^F(4)^1). 


In  this  standard  definition  we  have  also  used  a Stieltjes  integral,  mainly  to  abbreviate  the 
formula.  The  reader  who  is  not  familiar  with  these  will  have  no  trouble  if  he  remembers  that 
when  F is  differentiable. 


dF 

dF{4>)=^-d4>-, 


otherwise  dF(^')  is  the  concentrated  probability  » <^').  For  similar  brevity,  Stieltjes  integrals 
will  be  used  below,  when  this  is  convenient,  to  improve  the  display  of  the  relations  obtained. 

The  main  technical  advantage  of  Stieltjes  integrals  is  the  simplification  of  expectations 
involving  continuous  and  concentrated  pro^bilities.  When  the  latter  are  not  zero,  the  con- 
tinuous deiuity  will  be  termed  Improper,  retaining  the  same  notation.  Thus  in  general 

j*dF=i  but 

As  a practical  example  we  have  for  crack  lengths 

roo 

1 — I dF(alt)  = Probability  that  a structure  will  not  survive  to 
J 0 time  t (for  a crack  to  start). 

The  vertical  stroke  above  is  the  normal  convention  for  conditional  probability  so  that  a|t  reads 
as  “a  given  Without  attrition,  the  conditioning  event,  that  the  structure  hu  reached  a life  (, 
is  a certainty.  In  our  context,  however,  the  conditioning  event  is  this  or  prior  attrition;  /(a|()  is 
then  improper,  summing  only  to  the  reduced  probability  of  reaching  life  t. 


X COMHNID  CIACUNG  AND  ATTRIIION 

Consider  the  transition  of  possible  cracks  in  the  upper  part  of  Figure  3 over  a time  dr. 
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FIG.  2 RELATIVE  SCATTER  OF  INITIATION  AND  CRACKING  TIMES 
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FIG.  2 (cont)  CRACKING  IN  TWO  IDENTICAL  SPAR  BOOMS 


If  there  is  no  attrition  and/or /(a|0  it  the  density  of  crack  length  then  by  an  analogy  with  fluid 
flow*  the  constant  probability  that  cracks  will  lie  between  the  two  shown  may  be  expressed  in 
the  continuity  equation 

(3.1) 


where  £//Z>r  is  a total  derivative  along  the  crack  trajectory  and  for  single  cracks  dir  Ria)  >■ 
dRjda,  a known  function  of  a. 

Attrition  for  this  particular  case  will  correspond  to  the  risk  rate  r{a)  and  the  events  of 
attrition  and  cracking  are  mutually  exclusive.  If  the  deruity  / is  allowed  to  be  non-proper  then 
attrition  will  decrease  it  at  the  rate  /(a|r)r(a)  so  that  (3.1)  with  attrition  becomes 


with 


Dt 


(3.2) 


It  is  worth  noting  that  if  div  R is  retained  (3.2)  holds  for  several  cracks  but  r(a)  becomes  harder 
to  calculate.  Another  point  is  that  ria)  is  the  total  rate  of  attrition  including  thst  from  causes 
not  related  to  the  growth  of  fatigue  cracks  or  from  various  modes  of  failure.  Statistically,  the 
attrition  simply  adds  to  the  effect  of  the  dilatation  dR/da. 


FIG.  3 DERIVATION  OF  DEGENERATE  FOKKER- PLANCK  EQUATIONS 
AND  INITIAL  CONDITIONS 

3.1  nwrafterlarif  Eqaatleas 

Equation  (3.2)  is  a linear  first  order  partial  differential  equation  and  is  also  the  Fokker- 
Planck  equation  for  the  Markov  process  represented  by  the  successive  states  of  damage, 
crack  len^  and  total  attrition  as  a fonction  of  crack  length.  (The  usual  second  order  term  is 
abaent  because  (3.2)  is  based  on  deterministic  cracking  without  local  randomness.) 

From  well  known  results^*'^*'**  the  characteristic  equations  ate 

dt  - da/Ria)  - -4flf[dRlda  + r(a)]-  (3-3) 
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3.2  UHal  CoadMoas 

Equation  (3.5)  is  the  complete  formal  solution  of  (3.2)  but  the  initial  density /«  ^f{ao\to) 
is  still  to  be  found.  It  follows  from  the  same  type  of  continuity  considerations,  but  it  must  now 
be  remembered  that  the  density  of  initial  failures  is  reduced  in  general  by  attrition  of  uncracked 
structures  (e.g.  static  oveiload,  hijacking,  etc.).  The  attrition  of  cracked  structures  is  already 
accounted  for  by  (3.2)  where  r{a)  includes  any  risk  of  interest,  as  well  as  the  failure  of  fatigue 
weakened  structures. 

Let  F{t)  be  the  probability  distribution  of  initial  lives  in  the  absence  of  attrition  but  in  its 
presence  suppose  that  the  proportion  of  uncracked  survivors  is  reduced  to 

Hm  - m-  (3.6) 

Among  these  survivors  there  are  two  mutually  exclusive  forms  of  attrition,  namely, 

(а)  Failure  or  removal  of  the  uncracked  structure  with  a risk  rate  ro  =>  r(0),  and 

(б)  Transition  to  being  a cracked  structure  for  which  the  risk  rate  is  the  same  as  it  is  without 
attrition,  i.e./(0/(l  - F(/)). 

For  the  survivors  (3.6)  the  usual  reliability  equation  then  takes  the  form 

+ (3.7) 

in  which  only  H(l)  is  not  known.  If  H{0)  = 1 the  solution  is 

H(t)  - e-'<^  (3.8) 


It  is  now  possible  to  find  initial  conditions  for  (3.5)  from  the  preservation  of  probability 
between  AB  and  CD  in  Figure  4.  The  probability  of  initiation  betwera  A and  B is 

Fr(Survival) . Pr(Init>ation|Survival) 

Hm  - mm  - exp[-/(ryr/(l  - F(0)J}  - (3.9) 

to  the  first  order.  This  is  not  obviously  the  same  as  the  probability  of  a csack  between  C and  D 
because  some  attrition  may  occur  during  the  time  dt  with  the  risk  Hde+y  If  survival  of  this  is 
included,  like  the  risk  in  (3.9),  and  equated  to  the  probability  of  an  initial  crack  between  C 
and  D then  the  preservation  of  probability,  with  allowance  for  attrition,  is  expressible  as 

H/exp(-iiao+)dl) . dt  =/(n«|/  + dt)R{ao)dt. 

The  term  for  attrition  here  is  of  second  order  and  the  initial  condition  reduces  to 

/M)  = H{t)miR(a,)  (3.10) 

When  this,  with  (3.8),  is  substituted  into  (3.S)  we  finally  obtain 


where  the  initial  conditions  are  expressed  on  the  axis  of  time  or  cycles  and  is  the  initiation 
time  corresponding  to  the  crack  a\t.  The  neglected  second  order  term  in  (3.10)  corresponds  to 
simultaneous  initiation  and  hijacking. 


33  Distribotion  of  Final  Failures 

First  these  are  considered  without  inspections.  For  a given  crack  a\t  the  risk  of  failure  is 
Ka)  so  that  the  probability  of  fatigue  failure  during  ((,  t + dt)  is  riii)f{a\t)dt,  remembering 
that  fia\t)  here  has  been  defined  to  allow  for  previous  attritions  At  time  t these  previous 
attritions  (the  chance  that  the  structure  will  be  overloaded,  hijacked  or  removed  btfore  initial 
cracking),  continue  at  the  rate  r,fixp(—rot)(\  — m))-  When  this  is  added  to  the  expected  attrition 
of  cracked  structures  the  density  of  lives  to  failure  from  all  causes  appears  as 


r(a)f(a\t)da 


<Kt)  = roe-'m 

J ao  (3.12) 

which  includes  infinite  cracking  rates.  The  notation  ^ is  a convention  to  indicate  the  presence  of 
attrition.  The  occurrence  of  infinite  cracking  rates  at  finite  lengths  is  the  case  of  runaway  cracking 
and  is  discussed  below.  The  two  terms  in  (3.12)  form  improper  densities  for  failure  of  uncranked 
and  cracked  structures  respectively.  Now  consider  the  secoi^  term  and  substitute  for  the  density 
from  (3.11).  This  corresponds  to  the  attrition  of  cracked  structures  so  that 


^t  n Cracking) 


and  for  autonomous  cracks  the  integrals  may  be  simplified  by  a change  of  variable.  Retrace 
e by  tm,  the  time  that  the  crack  requires  to  grow  from  a«  to  c.  This  is  defined  for  an  autonomous 
crack  and  a moment’s  thought  will  convince  the  reader  that  da/dia  » R(a).  With  these  changes, 
and  the  hijack  term  of  (3.12),  (3.13)  becomes 


In  the  integral  the  factors  depending  on  r - /•  form  an  improper  density  of  initiation, 
since  the  first  term  allows  for  prior  attrition  by  ultimate  loads.  The  leinaining  factors  could 
obviously  be  the  result  of  attrition  with  the  risk  r(u)  ■ r<e(/a))  on  the  set  of  structures  which 
start  to  crack  at  ta  •*  0;  they  are  the  probability  density  of  the  time  spent  with  cracks  in  the 
structure  before  final  failure.  If  r«  « 0 then  the  convolution  integral  in  (3.14)  would  be  a proper 


Too  / 

\ 

1 - la)  exp  ( - 

litldf], 

Jo  \ 

Jo  / 

density  leflecting  the  fact  that  the  total  life  is  the  sum  of  the  initiation  period  and  an  independently 
distributed  time  of  crack  growth. 


34  Menwt  GwsraHng  Fnetton  for  Lives 
With  (3.14)  the  life  distribution  (3.12)  becomes 

i (/)  = r,e-'<*  (1  -FiO)  + e-'o^fO,)  •/(/.) 
so  that  the  Moment  Generating  Function  Eexp(ut)  is 

Af#(«)  = (1  - Fit))dt  + /(j  - t,)dF{U 

= ^(1  - F(f ))  de<*  j*  J ~ “** 


= Mdu  - r.)A/.(u)  + — —(Mdu  - r.)  - 1)  (3.16) 

U — To 

where  Mt  and  Mt  refer  to  the  initiation  and  cracking  (imes  respectively.  When  there  are  no 
ultimate  loads  ro  0 and  the  independence  of  t and  u is  confirmed.  The  result  is  exact  for  autono- 
mous cracks  and  some  expansions  of  the  cumulants  of  final  life  ate  shown  in  Table  1. 

3J  Allowance  for  Inspections 

The  effect  of  an  inspection  may  be  regarded  as  a concentrated  attrition  applied  instan- 
taneously. In  terms  of  the  Fokker-Planck  Equation  (3.2)  and  Figure  5 the  boundary 
conditions  are  altered  for  different  crack  lengths  along  the  inspection  boundaries  t Tt.  The 
computation  of  /(a|r)  and  therefore  of  4(t)  pro(%eds  in  turn  from  regions  /i  to  /m 

with  changes  of  boundary  conditions  at  Tt. 


Inspection  boundaries 


Initiation 

boundary 


I \ \ Domain  of 

Domain  of  depandenoa 

dependence  on  AB 

on  f (t),  T,  < t < T,  ^ 


Number  of  cycles  or  time 
FIG.  5 DEPENDENCE  ON  INSPECTIONS 

Let  P{a)  be  the  operating  characteristic  of  the  method  of  inspection  so  that 
Fifl)  Pr  (Not  acting  about  a crack|Laagth  «). 

In  general  P(«)  will  be  the  result  of  operating  experience  but  the  movs  exotic  possibilitiss  of 


finding  cracks  and  not  acting  or  randomly  rejecting  structures  can  be  included  in  subsidiary 
calculations  for  it.  After  the  first  inspection  the  boundary  conditions  along  h become 

/(airi+)  - mmn-) 

but  this  will  only  affect  the  cracks  that  have  already  begun,  that  is  those  in  the  domain  of 
dependence 

(a,  0 e h and  a ^ a(/  - Ti) 

which  is  that  part  of  A above  the  characteristic  crack  trajectory  starting  at  (a*,  Ti)  (Fig.  5). 
Similarly  the  trajectories  from  the  other  (a«,  7|);s  will  divide  all  the  It  into  other  domains  of 
dependence  which  depend  on  conditions  set  at  the  inspection  boundaries  t Tt,  a ^ a,  or  else 
the  initiation  boundaries  Tt  <t  ^ Tm,  a ~ a,.  Along  each  of  the  latter  the  probability  of 
initiation  will  be  greatly  affected  by  the  inspection  policy.  The  study  of  these  policies  and  their 
effect  is  related  to  renewal  theory.^*  Because  little  is  known  at  present  and  the  subject  promises 
to  be  vast  it  will  not  be  pursued  here  except  for  the  mention  of  four  typical  possibilities.  These 
are  o 

(a)  the  structure  is  removed  from  service; 

(b)  it  is  repaired  as  new,  i.e.  the  damage  is  the  same  as  that  of  a new  structure; 
and  finally  there  are  the  more  difficult 

(c)  repaired  to  a = 0 but  with  residual  fatigue  damage  (this  means  that  the  repair  will  tend 
to  fail  more  rapidly  than  originally); 

(d)  replaced  by  new  structure  selected  at  random  from  the  population. 


4.  EXTENSIONS 

For  the  model  above  the  obvious  need  is  more  investigation  of  inspection  policies.  However, 
there  are  other  simple  extensions  which  are  either  interesting  or  related  to  current  practice. 

The  first  of  these  is  that  the  initial  crack  length  a«  could  be  taken  as  zero  in  the  context  of 
the  present  development  although  there  are  physical  reasons^'*’*  why  it  should  not  be. 

4.1  KaMway  Cracks 

In  Forman,  Kearney  and  Engle's  relation^*  for  daldi  and  for  practical  purposes  in  some 
experimental  results,  it  is  possible  for  R(a)  to  become  infinite  at  some  finite  length  c,  after  the 
. growth  time  A say.  If  this  happens  (defining  a “runaway  crack”)  then  obviously  such  structures 
survive  no  longer;  in  the  basic  statistical  terms  r(c)  1 and  it  will  transpire  that  e -*co  makes 
no  difference.  The  effect  of  this  is  to  telescope  all  of  the  failures  which  would  have  occurred 
later  on  (had  daldt  not  become  infinite)  into  ffie  instant  at  u + A- 

This  is  the  continuous  idealisation  of  the  final  cycle  in  a discrete  sequence  of  loads. 

4.1.1  Discrete  Analogy 

It  may  be  argued  that  a risk  approaching  r(c)  when  t < t»  + te  would  leave  no  survivors 
for  such  a wholesale  demise  when  the  risk  become  unity.  However,  we  shall  demonstrate  the 
effect  for  (diysically  more  realistic  discrete  time,  i.e.  individual  load  cycles.  It  is  suflkient  to 
postulate  a fixed  crack  always  “beginning”  at  the  origin  with  the  induced  risks  r»,  n . . . rH_i 
and  of  course  TN  = I. 

Then  forj  > N ~ 1 

i>r(Survivingy  cycles)  - (1  - r,) ...  (I  - n)  - (1  - Ff),  say  (4.1) 

In  particular  the  nature  of  r/  ensures  strictly  that 

0 < 1 - F»_i<  1.  (4.2) 

since  S is  finite.  The  attrition  on  the  next  cycle  is  complete, 

riKl  - Fii_i)  « 1 - F»_i, 
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and  the  range  of  the  life  distribution  is  (0,  N)  with 

/n  - 1 - /■»-!, 

the  telescoping  effect. 

If  this  simplified  case  is  repeated  in  continuous  time  then  (4.1)  is  replaced  by 


1 - F(tc)  = 


= exp  I “ 


as  in  (4.2)  if  the  integral  converges,  as  it  may. 

As  in  (4.3)  /(/<•)  is  a concentrated  probability,  an  exceptional  case  in  continuous  time,  but 
only  because  (4.4)  is  an  approximation  to  (4.1).  It  is  easily  shown  that  there  exists  a believable 
r(0  for  which  (4.4)  uniformly  exceeds  (4.1). 


Runawmy 
limit  of-> 


May  crack 
of  A.  B.  / 

I rx/ 


Critical  crack  length 


Cyclat  or  tima 


FIG.  6 RUNAWAY  CRACKS 
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In  this  section,  runaway  cracks  will  be  appraached  as  the  limiting  case  of  an  inequality 
for  the  distributioa  fiinction  of  crack  lengths  F(a|l).  For  this  purpose  suppose  that  c ■■  4r«)  is 
some  arbitrary  length  above  which  the  rate  is  so  that  f(ao)  - t{e)  is  “short”. 

In  the  notation  of  (3.IS)  let  ns  divide  the  infinite  range  in  (3.14)  at  die  arbitrary  times  r« 
and  Ta  > r«.  The  first  of  theee  will  become  the  time  until  a ciack  runs  sway  whibt  7#  is  to  allow 
for  the  infinite  limit  in  (3.14). 


Thus  (3.14)  may  be  written 


^tnO 


- O{e-roUf(0^\  _ F(r.))] 

where  0 > / > T — tc  by  the  mean  value  theorem.  Because  /(r)  is  zero  when  t <0  this  inequality 
can  be  sharpened  to  0 ^ f < f — te  since  Pr{i  + u > l)  is  then  zero. 

The  third  term  is  essentially  positive  and  its  neglect  indicates  that  always 


with  equality  when  fCT,)  - 1.  This  happens  when  Ta  is  infinite  or  at  finite  values  when  r(c)  » 
1,  i.e.  there  is  runaway  cracking.  This  is  the  case  we  discuss.  In  (4.5) 

M-r,f+  In(/(r  -tc-  O/fO  -tc)} 

- rof  + ln{l  - if ’If)  ^ cm-  i(r,  -f'/f(t  - h)). 

So  far  the  introduction  of  f or  Af  has  avoided  the  use  of  unit  risk.  In  this  Section  risk  must 
be  correlated  with  cracking  so  that  ric)  -*■  1 is  implied  by  f,Af  -*■  0,  i.e.  certain  failure  follows 
instantaneous  fracture. 

For  runaway  cracks  in  continuous  time  then  n C)  becomes 


From  Section  3.3 


and  in  the  limit  (3.12)  or  (3.1S)  take  the  most  general  form 


The  last  term  is  the  “risk  of  fatigue  fracture”  as  defined  by  Payne  and  Diamond.*  The  moment 
generating  function  (3.16)  is  not  altered  but  if  X(a(te))  is  infinite  then  Afa(u)  needs  to  be  computed 
using  the  concentrated  probability  (4.6)  at  /«. 


41.3  Msrrstsneas  Eirett 

In  (4.5)  txpM  - I + M indicates  the  effect  of  fracture  not  being  instantaneous.  If  r(c)  > 1 
then  f cannot  “exceed”  one  load  cycle  and  with  this  extreme  bound  we  now  show  that  M is 
negligible. 

Suppose  that  unit  time  allows  an  average  of  in  cydee  so  that  heuristkally 
r.-mr«(0) 

where  the  subecript  m indicates  discrete  load  quantities.  Substitution  with  I ~ 1/m  than 
indicates  that  M < rij(0)  - - r«)(l  - F(t  - f«)]  which  is  small  for  practical  casas. 


Too  / 

ftc  \ 

dF(u)  = exp  ( - 

- r(ta)dt,\ 

J ft  \ 

Jo  / 
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So  r«r,  the  crack  rate  has  depended  only  on  the  current  length  so  that  any  possible  growth 
curve  is  a translation  in  time  or  cycles  of  one  that  begins  at  the  origin.  In  practice,  the  rate,  say 
daldt  ■■  R(.a,  u),  depends  randomly  on  the  particular  specimen,  here  itemised  as  u.  For  a two 
model  of  fatigue  this  variability  can  either  be  neglected  in  comparison  to  the  scatter  from 
the  H«m«ge  phase  or  else  (4.7)  can  be  repeated  for  different  cases  and  averaged.  Within  the  present 
model  the  crack  rate  can  also  be  modiM  to  allow  approximately  for  random  crack  rates. 

In  other  cases,  notably  for  some  welded  structures,  it  may  be  more  realistic  to  auume  that 
a crack  begins  immediately  or  at  some  fixed  time  and  another  further  idealisation  would  be  that 
the  time  to  a given  crack  length  has  some  given  probability  density .<•* 

Such  cases  are  included  by  taking  a crack  rate  Jl(a,  /)  depending  both  on  time  and  length 
which  also  generalises  the  two  sUge  model.  Equation  (3.2)  still  holds  together  with  the  subsequent 
argument  except  that  the  dilatation  of  R is  now  a partial  derivative  depending  on  time.  The  first 
effect  of  this  is  that  (3.S)  fails,  which  in  turn  means  that  the  simple  form  (3.16)  for  the  moment 
generating  function  must  be  abandoned;  the  propagation  time  and  initiation  are  not  independent 
However,  some  salvage  is  possible.  Crack  growth  trajectories  are  still  characteristics  and 
along  these 

d.  , iR(a,  t)[ia  -f  r(d) 


leading  to  a formal  solution  with  the  previous  boundary  conditions.  Consider  the  total  derivative 
of  R along  the  crack  trajectory. 

Da  da  it  da . 

Replacing  iRfia  in  (4.8)  and  integrating 


Initially 


t»  depending  on  a. 


which  corresponds  to  (3.1 1).  The  integrals  here  are  defined  by  relating  a and  t along  the  charac- 
teristics; (3.12)  follows  as  before  with  the  changes  above. 

The  modifications  needed  to  display  runaway  cracking  may  also  be  extended  to  non- 
autonomous  cracks.  The  argument  is  similar  but  must  begin  at  (3.12)  which  now  becomes 


Since  R(e)  is  infinite  the  last  term  becomes  unity  for  the  range  of  integration  (c,  eo)and  thie  leavee 

^ -/waMi 


as  a probability  density  from  an  attrition  along  the  characteristic  through  (c,  t)  with  the  risk 
rate  r/X.  In  addition  there  is  no  longer  any  reason  why  e need  be  independent  of  r.  All  of  these 
considerations  lead  to 


/(t^a))  exp 


+ - c)), 

as  the  most  general  one-crack  form  of  (4.2),  where 


along  the  crack  tnuectory  through  (a,  i)- 


43  Diati«irtiaaarFaHi«Loaia 

In  the  context  of  this  paper  failing  loads  may  be  inferred  from  the  crack  length  when  the 
structure  fails  so  that  we  address  the  problem  of  finding  the  density  of  crack  lengths  at  failure. 
With  respect  to  the  overall  probability 

Ku)/(a|0  "■  PKAttrition  n Survival  to  t n Crack  length  a). 

By  definition 

4(0  **  ?r(Attrition  at  r) 

in  overall  measure  whence  the  crack  length  density  conditional  upon  attrition, 

4(ol0  - r(a)/(alOm  (4.11) 

at  a given  number  of  cycles; /(a|0  refers  to  cracks  which  may  be  merely  present  at  time  /.  For 
autonomous  cracks  and  > a > c,  substitution  from  (3.1 1)  produces 


4(tmi0  - e-roto/(,,) 


« e-f.'‘-w/(r  - (4.12) 

with  the  change  of  variable  a-*!,  used  for  (3.14). 

In  words,  the  bivariate  density  of  life  and  crack  length,  for  life  at  final  failure,  is  the  product 
of  growth  time  density  and  the  density  of  initiation  of  the  characteristic  divided  by  the  current 
rate  of  cracking. 


4J.1  ZareCkucha 

The  density  (4.12)  should  be  augmented  by  concentrated  probabilities  that  the  cracks  in 
the  failed  structure  eitlier  did  not  exist  or  possibly  were  runaway  cracks. 

The  first  of  these  is  the  conditional  probability  ^0|0  of  atbition  by  the  risk  r*.  From  the 
first  term  of  (3.13)  or  (4.I0X  normalised  by 

^r)«(0|0  - - ^0)  (4.13) 

wUdi  may  be  integrated  u it  stands  for  dis  marginal  probability  of  fiulure  or  removal  widiout 
craddng.  Partial  integration  with  respect  to  the  Actor  in  r*  leads  to 


(4.14) 


♦(0)  - 1 - MK-r,) 
- rtui  + 0(r,*), 


where  m - Mean  cycles  to  initiation. 

It  is  suggested  that  (4.14),  the  product  of  arithmetic  mean  life  to  initiation  with  risk  rate  of 
uncracked  structures,  is  an  appropriate  measure  of  "fatigue  sensitivity”  since  it  is  the  probability 
that  failure  will  not  be  related  to  fatigue. 


4J.2  RBMway  Oacks 

By  arguments  the  same  as  above  the  last  term  of  (4.7)  may  be  normalised  to  the  conditional 
probability  of  runaway  cracking  whence 

^/Xl  - <l»(c|0)  - e-«‘ -*'*/(<  - /.XI  - - c))  (4.15) 

with  /(/•)  from  (4.6).  Since  /(/)  is  zero  for  negative  times  this  integrates  over  f to  the  marginal 
prolMbility  of  runaway  cracking 

(1  - dKc))  - MX-r.Xl  - f=(/.|a  - c)) 

(4.16) 

« (!  - r.,i,)(l  - HU\a  - c)). 

When  me)  is  not  necessarily  infinite  but  large  enough  to  be  considered  dangerous  (4.16) 
may  be  regarded  as  a rough  measure  of  “fatigue  danger”.  The  interpretation  is  precise  for  run- 
away cracks  but  otherwise  clouded  by  considerations  such  as  those  of  Section  4.1. 
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When  (4.12)  is  integrated  with  respect  to  / - /«  one  obtains  the  nmrginal  density 

^a)  - WX-r,)/(/.)/«(o)  (4.17) 

which  follows  from  (3.11)  and  (4.11).  Inspection  of  these  original  equations  indicates  that,  just 
as  /(/•)  is  the  outcome  of  an  attrition  process  with  the  risk  r(a)  and  a a(U),f(le)IItio)  is  both 
the  transformation  of  that  density  to  that  of  a(U)  and  the  result  of  attrition  along  crack  length 
with  the  risk  r(o)//l(a)  (see  Section  3.3). 

It  can  be  shown  that  (4.14),  (4.16)  and(4.l7)provide  the  complete  marginal  distribution  of 
length  at  failure. 

5.  mSCUSSiON 

The  preceding  evolution  of  two  stage  fatigue  depends  only  on  the  density  of  initial  life 
f(t)  and  the  average  crack  rate  Jt(a).  If  failure  is  included  a risk  r{a),  averaged  over  streu  sequence, 
material  toughnen  and  untoward  circumstance  is  also  needed.  These  three  flinctions  then 
completely  define  fatigue  and  life  of  the  structure.  They  may  be  estimated  by  any  standard 
method,  experimental  or  theoretical. 

The  computer  programs  may  also  allow  for  inspections  if  the  operating  characteristic  ^n) 
of  the  overall  inspection  method  is  known. 

Most  current  experimentation  about  crack  rates  tacitly  assumes  that,  apart  from  the  stress, 
they  depend  only  on  the  configuration  of  the  crack  at  the  time.  This  is  the  aaumption  of  autono- 
mous cracking  which  is  suited  to  the  use  of  fracture  mechanics  and  support^  by  the  small 
selection  of  data  in  Figure  2 and  most  of  that  in  the  literature. 

For  single  cracks  we  have  shown  that  the  assumption  of  autonomous  cracking,  which  is 
essentially  that  fatigue  has  two  stages,  leads  to  remarkable  simplificationt  of  the  associated 
rsUaWlity  theory  provided  that  tUs  is  incorporated  into  the  structural  IhtigHe  analysis.  There  is 
a parallel  increase  in  understanding  of  the  interaction  of  ultimate  loads  and  the  fhtigue  of  stme- 
tarse  which  may  be  summarised  by  Equations  (3.16)  and  (4.7). 

Both  of  these  are  exact  for  a very  general  model  with  automatic  aRowance  tat  the  Increasing 
strength  of  structures  with  time  as  the  weaker  onm  are  weeded  out  by  higher  risks.*  Some  iatar- 
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mation  about  this  will  be  contained  in  the  crack  length  distribution  results  of  (4.12),  (4.13)  and 
(4.  IS)  though  interpretation  in  terms  of  strength  depends  on  the  nature  of  the  subsidiary  problem 
of  estimating  the  risk  lia).  In  the  same  section  it  can  be  seen  that  the  marginal  distribution  of 
strength  particularly  the  virgin  strength  and  that  for  runaway  cracks  provides  natural  measures 
of  the  sensitivity  of  the  structure  both  to  fatigue  in  general  and  to  fatigue  which  may  be  regarded 
as  dangerous. 

Finally,  many  of  the  results  can  be  extended  to  cracks  that  are  not  autonomous  with  no 
lou  of  accuracy  and  very  little  loss  of  simplicity. 

The  representation  of  the  crack  length  density  as  a Markov  process  governed  by  a partial 
differential  equation  also  allows  an  exact  solution  for  the  results  of  general  inspections  with 
realistic  operating  characteristics.  Although  this  is  exact  in  principle  the  calculations  involve 
following  boundary  conditions  across  many  contiguous  domains  of  dependence.  For  single 
cracks  the  optimisation  of  such  computations  for  either  speed  or  accuracy  is  the  major  question 
in  the  application  to  structural  fatigue  and  reliability. 

For  structures  with  several  cracks  the  success  of  the  results  here  demonstrates  that  any 
theory  of  structural  fatigue  must  intimately  include  reliability  and  the  latter  subject  must  now 
become  multi-dimensional.  In  practice  many  structures  fail  by  a single  dominant  crack  and,  if 
this  is  so,  many  of  the  present  results  will  apply  if  the  initial  crack  length  a,  is  taken  to  be  some 
length  exceeding  that  of  non-dominant  cracks. 

When  they  coalesce,  non-dominant  cracks  will  thus  affect  the  stress  distribution  and  qualify 
for  our  definition  of  a crack.  Before  this,  however,  their  progress  would  be  described  as  con- 
tributing to  damage. 
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